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1. Introduction 

A very useful method of studying linear operators by decomposing the space on 
which they act into invariant subspaces is known as spectral theory. An example 
of an application of spectral theory is the problem of finding a set of eigenvectors 
or diagonalizing a linear map on an infinite-dimensional space [HN]. When a finite 
dimensional linear operator is diagonalized, there exist a set of eigenvalues and their 
corresponding eigenvectors. Along the directions of an eigenvector with its given 
eigenvalue, the action of the operator is just multiplication by the eigenvalue. The 
spectrum contains the set of eigenvalues that is also called point spectrum. In infinite 
dimensional case, the structure of the spectrum will often be more complicated such 
that there may exist a continuous spectrum or residual spectrum which do not contain 
a set of eigenvalues and are different from the point spectrum. In my thesis, I study 
the spectrum of an operator based on some numerical results of it and also properties 
of the discrete Laplacian operator. 

The discrete Laplacian operator A is defined on £^(Z) by 

{Ax)k = {xk-i + Xk+i) - 2xk, 

where x = {xk)'kL-oo ^ ^^(Z,). Moreover, three features of A are described in a 
theorem as follows: 

Theorem 1. Let S be the right shift operator and S* be the left shift operator. A = 
S + S* — 21. The spectrum of A is entirely continuous and consists of the interval 
[-4,0]. 

Proof. First, ||A|| < ||5|| + ||»S*|| + 2||J|| = 4 where || • || denotes the operator 
norm. Second, A = —{S — X){S*—X) = —{S — X){S—X)* < 0. These two arguments 
imply that cr(A) C [—4, 0] where cr(A) denotes the spectrum of A. The facts that 
(t(A) = [—4, 0] and (7(A) is purely continuous will be shown in Section 3. 

The discrete Laplacian operator is one of the most important and oldest difference 
operator [D] which is closely related to the operator I'm studying in this paper. But 
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before I introduce this operator, I would like to talk about Schrodinger's equation 
and Schrodinger operators. 

A Schrodinger's equation without a potential term is a partial differential equation 
defined as = —Au [E]. On the other hand, Schrodinger operators acting on L^(M) 
are operators such as Hf{x) — —Af{x) +v{x)f{x), where F is a real- valued function 
on R and is called a potential [D]. We may impose condition on V such as choosing 
V E L^(M) so that H is an self-adjoint operator on L^(]R). For a detailed proof of 
this specific condition, see Davies [D]. Moreover, '/ e L^(M) with ||/||2 = 1 is called 
a wave packet or state, and represents the instantaneous configuration of a collection 
of electrons, atoms and molecules. The operator H is also called the Hamiltonian for 
historical reasons - quantum theory can be regraded as a non-commutative version 
of classical Hamiltonian mechanics. The evolution of a quantum system is controlled 
by the Schrodinger's equation — Hf with solution f{x,t) — e'^^^ f{x,Oy [D]. 

Finally, I want to now introduce the problem of my project. The term 'A' which 
you will see below is a coefficient, not the discrete Laplacian operator mentioned 
above. 

For a particular example of a Schrodinger's equation which is also the main subject 
of my thesis, a time-dependent discrete Schrodinger's equation defined on £^(Z) can 
be written as follows: 



^ cosh{r]{x — 1 — r)) cosh(77(x + 1 — r)) 
with A = cosh(77), 77 e R+, x e Z, r e R. 

In order to understand the operator i^, we define and study H : £^(Z) — > £^(Z) as 
follows 



d 
I— 
dt 




X 



where 



2 cosh{rj{x — r))^ 
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such that 
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a, is hnear and bounded. The fact that V, is hnear is trivial. Ti is bounded because 
ll^ll < + 11^5*11 + ||e.J|| < ||5|| + ||5*|| + ||2J|| = 1 + 1 + 2 = 4 since the maxima 

of 1/A and are 1 and 2, respectively. The reason that the maximum of is 2 for 
any r and 77 is because of the properties of which will explained in Section 2, and 
also Figure 1 below plots an example of versus x. 

As you can see that the two main differences between the discrete Laplacian op- 
erator and H. are the tridiagonal terms of their matrices where those terms in the 
discrete Laplacian operator are constant functions: 1, —2, and those in Ti are vari- 
ables depending on r and 77: —1/A, respectively. 

We know that the spectrum of the discrete Laplacian operator contains only the 
continuous spectrum when the diagonal terms of the operator are constants —2. On 
the other side, the diagonal terms of the operator Ti. have a trapping region around 
their centers and are 2 for large x. Figure 1 presented below shows an example of 
with certain parameters r and rj. So, 7i is a perturbed case of the discrete Laplacian 
operator, and the nature of the structure of the spectrum of H. is different from it of 
the discrete Laplacian operator because of the term e^- 

In this paper, I study the spectrum of H. defined above based on the discrete Lapla- 
cian operator which is an unperturbed case of H. and also some numerical results that 
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Figure 1. e^^ vs x with r = 1 and rj = 1 

can be used to investigate the solutions of the time- dependent discrete Schrodinger's 
equation from the following three aspects: 

• stationary solution (Tiz/ = 0) 

• periodic solution {Hu = \u) 

• properties or representatives of the general solution 

I have only studied the stationary and periodic solutions but not yet the properties 
or representatives of general solution here. 

For infinite dimensional Ti depending on r and r], my research shows that there exist 
three eigenvalues depending on a certain region of r and rj. The first eigenvalue of Ti 
which equals zero is exact. A detailed analysis regrading this eigenvalue also called 
the zero mode [MN] will be discussed in Section 3. Moreover, a second eigenvalue 
exists which is proved by Michoel and Nachtergaele [MN]. Finally, we also predict 
that a third eigenvalue exists for certain parameters r and r] based my numerical 
results. The primary focus in my thesis is to find critical values or intervals of the 
parameters r and rj where these eigenvalues exist by numerics. More precisely, the 
main object is to determine the particular regions where there exist one, two, or 
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three eigenvalues for infinite dimensional Ti.. In order to study this operator by 
numerics, we concentrate on a finite dimensional subspace of ^^(Z) such that is a 
finite dimensional operator. Since there always exists a set of eigenvalues for finite 
dimensional we therefore focus on the eigenvalues that are isolated from the set 
of all other eigenvalues which will belong to the continuous spectrum and consider 
those isolated ones as the candidates of eigenvalues which will belong to the point 
spectrum for infinite dimensional Ti. Those 'isolated eigenvalues' I mentioned can be 
seen more clearly on Figure 3 and Figure 4 in Section 3. 

My numerical results are created by three programs written by myself in MATLAB. 
First, I graph the spectrum of Ti. versus both parameters r and 77 separately to see how 
the set of the point spectrum is distributed with different r and 77. Second, I graph the 
eigenvectors with respect to the three eigenvalues versus their components to see how 
the eigenvectors are distributed and I'm hoping to see that the eigenvectors are still 
nonzero for large matrix size n with certain r and 77. Third, I calculate the maxima 
of the square normalized components of all those three eigenvectors. I compare the 
maxima with the same parameters and difi^crcnt matrix sizes individually and my 
expectation is that if the maxima do not decrease at certain values or intervals of r and 
?7, then the eigenvectors are not zero which also implies that there exist eigenvalues for 
infinite dimensional 7Y. This step mainly support my results. Finally, for the third 
eigenvalue mentioned above, I also observe three components of its corresponding 
eigenvector to see the decreasing rates of them. My meaning of decreasing rates 
will be explained in the next few sections. This process is to determine whether my 
conclusion of the existence of the third eigenvalue at certain r and rj based on the 
previous steps are correct or not. The next three sections show my analysis in detail 
of the problem and the last section concludes my entire paper, and they are presented 
in the order as follows: 

Primciry Results: This section states the procedures and the main discoveries 
in my project. 
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The Spectrum: This section discusses the definition of spectrum and also the 
structure and feature of the spectrum of H.. 

The Eigenvalues: This section shows my numerical results for all the eigen- 
values that are found. 

Conclusion: This section summarizes my results and also talks about some 
further discussions of this operator. 
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2. Preliminary Results 



If 7i is a finite dimensional operator, then it is also a compact operator. A compact 
operator is bounded, and a compact operator that is also symmetric is self-adjoint. 
For more details, see Hunter and Nachtergaele [HN]. Therefore, H is a, self-adjoint 
operator. Moreover, ?i is a nonnegative operator. In order to prove this fact, we need 
to use the following lemma first. 

Lemma 1. Let A be a 2 x 2 matrix such that 



where a,b & M., c & C Then A>0 if and only if a,b >0 and ab > |cp. 

Proof. Let ^ > 0. First, ^ > if and only if all its eigenvalues Ai,A2 > 0. 
Second, the determinant of A after simplification is — {a + b)X + ab — |cp = 0. 
Moreover, 



Since Ai,A2 > 0, we have a6 > |cp > and a, 6 > 0. Conversely, if a, 6 > and 
ab > |cp, then Ai, A2 > which also implies that ^ > 0. 

Now, we can decompose 7i into infinitely many submatrices such that 7i = X^^gg hx 
where 





A1A2 — ab — \c 



















1/A 



1/A e, 
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cosh.{ri{x — r)) 



A cosh(?7(a; ± 1 — r)) 
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Since only four entries in are nonzero, if we want to prove that is nonnegative, 
it is sufficient to focus on a 2 x 2 matrix hx defined as follows: 

\ -i/A J 

Since > for all x e Z, > and €~_^.l > 0. Thus, the remaining is to show 
that whether e^e^_^_i > (1/A)^ is true or not. The calculation is presented below. 



cosh(?7(,T — r)) 
' A cosh(?7(a; + 1 — r)) ' 



cosh(r/(a; + 1 -r)) _ 3 
^ Acosh(77(x-r)) ' ^ ' ' 

= 

Hence, hx is nonnegative for all a; G Z because e+,e~^^ > and e1^e~^i = (1/A)^. 
So, hx is nonnegative for all x & Z which implies that H is nonnegative and also all 
eigenvalues of Ti. are nonnegative. 

For infinite dimensional Ti, my numerical results predict that there exist three iso- 
lated eigenvalues of Ti. for a certain range of the parameters r and rj such that the point 
spectrum of H is not empty, unlike the discrete Laplacian operator. My conclusions 
are mainly based on calculating the maxima of the square normalized components 
of all the three eigenvectors corresponding to their eigenvalues with various r and rj 
and also observe the pictures of all the three eigenvectors and the decreasing rates for 
three chosen elements of the eigenvector of the third eigenvalue to compare them with 
my numerics because the range of the existence of the third isolated eigenvalue with 
respect to r and r] is different from it of the other two isolated eigenvalues. Roughly 
speaking, with certain interval of 1], the first two eigcnvahics exist for all r in a certain 
region which will be determined but the third eigenvalue does not exist for all r in 
that region. More details are presented in the following few paragraphs. 

I am using finite dimensional H to approximate infinite dimensional Ti. The rea- 
son why I follow such path is based on what we know about the discrete Laplacian 
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operator. If 77 = 0, then Ti is different from the discrete Laplacian operator only by 
a sign. Moreover, my numerical data includes the case where 77 = and it shows 
that the eigenvectors approach to zero when the size of Ti becomes larger and larger 
though not presented here. So, we expect that all the eigenvectors of H. when 77 = 
converge to zero by the sup norm as the size of Ti. goes to infinity which indicates that 
the hmits of the eigenvalues are in fact not eigenvalues. The reason that the converge 
is in the sense of sup norm is because my numerical data shows that the maxima of 
the square normalized components converge to zero, and so the convergence is also 
in the weak £^ sense. Moreover, we already know that the spectrum of the discrete 
Laplacian operator does not contain any eigenvalue. Therefore, using finite cases of 
H to predict infinite cases of Ti. is reliable even if we change the parameters r and r]. 

Now I want to talk about why I start with certain r and 77 out of other possibilities. 
For an n-dimensional Ti where ti e N, the diagonal terms (ej;)"^_„ are exponentially 
localized functions, depends on r and 77 in which r plays the role of shifting and 
77 concentrates ex around its center, and Cx is also symmetric around its center shown 
in Figure 1 above. It is therefore sufficient to focus on the case r G [n/2,n/2 + 1] 
where n represents the size of Ti. Moreover, Ex is an even function because it is 
governed by hyperbolic cosine functions which are even. As a result, we can begin 
with 7/ = and then gradually increase rj. By plotting the graphs for the set of 
eigenvalues of Ti versus the parameter r G [n/2,n/2 + 1] with increasing values of 7] 
which will be displayed in the next section, the range of the set of the spectrum of 
Ti becomes smaller and smaller and three isolated eigenvalues are extracted from the 
spectrum. For infinite dimensional case, these three eigenvalues will belong to the 
point spectrum and all others will belong to the continuous spectrum. The concept 
and structure of spectrum will also be discussed in the following section. 

From my numerical results, there always exist two eigenvalues for all r G [7i/2, 7i/2+ 
1] and 7/ > 0. Moreover, for (r - [r]) G [0, 0.1] U [0.9, 1] and 7/ > 1, and for (r - [r]) G 
[0.2, 0.3] U [0.7, 0.8] and 77 > 2, there exist a total of three eigenvalues. The graph 
shown below predicts the general results for the locations of the existence of the 
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three isolated eigenvalues with different r and 77 for infinite approximation of Ti. For 
numerical data which is not completely presented here, the third eigenvalue exists 
when (r — [r]) 0, 1 and = 1; (r — [r]) 0.15, 0.85 and r] = 1.5; (r — [r]) ^ 0.25, 0.75 
and 1] = 2, and the two curves in the graph are produced based on connecting those 
points. More numerical results are therefore needed to verify this graph. 




0.5 - 

qI 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Figure 2. rj vs [r — [r]), where [r] denotes the integral part of r 

The results above show that the point spectrum of infinite dimensional Ti is nonempty 
at these certain regions. Hunter and Nachtergaele [HN] and Michoel and Nachtergaele 
[MN] have proved the existence of the first and second isolated eigenvalues analyti- 
cally, respectively. My results mainly come from comparing the maxima of the square 
components of normalized eigenvectors with different matrix size n at certain values 
of r and 77 numerically. The idea is that when the maxima do not decrease for certain 
r and f] while the matrix size n becomes larger and larger, the eigenvalue does not 
vanish because its corresponding eigenvector is nonzero. More details of my numerical 
calculations of the results presented above will be discussed in the next two sections. 

My numerical results including the maxima and the graphs of eigenvectors are 
created by a program written by myself in MATLAB named "mainresults.m" and 
the one for observing the decreasing rates is called "decreasingrates.m" . Also, the 
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pictures of the spectrum of H. versus (r — [r]) and 1/A, respectively, are generated by 
" spectrum, m" . The purposes of these three programs will be explained in Section 3 

and Section 4. The next section explains the spectrum of 7i in more details and also 
discusses the spectrum of the discrete Laplacian operator. 
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3. The Spectrum 

In some circumstances, there are no eigenvalues for an infinite dimensional bounded 
linear operator and so it is not possible to expect to find an orthonormal basis that 
consists entirely of eigenvectors. Thus, we need to define the spectrum in a more 
general way, instead of considering it only contains eigenvalues [HN]. 

The following two definitions are rewritten from 'Applied Analysis' by Hunter and 
Nachtergaele [HN]. 

Definition 1. Let A be a bounded operator defined on an infinite dimensional Hilbert 
space. The resolvent set of A, denoted by p{A), is the set of complex numbers A such 

that {A — AX) is one-to-one and onto. The spectrum of A, denoted by o^{A), is the 
complement of the resolvent set in C, meaning that (t{A) = 'C\p[A). 

According to the open mapping theorem, {A — AX)~^ is bounded if ^ — AX is one- 
to-to and onto, and therefore ^ — AX and {A — AX)~^ are one-to-one, onto, bounded 
hnear operators when A G p{A). The following definition gives the structure of the 
spectrum of a bounded linear operator. 

Definition 2. Suppose that A is a bounded linear operator on a Hilbert space. 

(1) The point spectrum of A consists of all A e cr(^) such that A — AX is not 
one-to-one. In this case A is called an eigenvalue of A 

(2) The continuous spectrum of A consists of all A e '^{A) such that A — AX is 
one-to-one but not onto, and the range of A — XI is dense in this Hilbert space. 

(3) The residual spectrum of A consists of all A e cf{,A) such that A — AX is 
one-to-one but not onto, and the range of A — XI is not dense in this Hilbert 
space. 

We define 7^ on a finite subspace of £^(Z) for our problem and we then have T-Lv — Xu 
where A is an eigenvalue and u is a nonzero eigenvector. Thus, the kernel of H does 
not contain zero, and so H — AX is not one-to-one and A e a{7i). 
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Because of the diagonal term in 7i, the structure of the spectrum of Ti is different 
from it of the discrete Laplacian operator. As shown in Figure 1 above, the function 
ex remains constant for large |x| but drops toward zero for small \x\ around its center. 
The existence of this dropping region makes the spectrum of Ti so special in which 
the point spectrum is not empty such that we expect there exist eigenvalues when Ti 
is infinite dimensional. The location, width, and depth of this trapping region depend 
significantly on the two parameters r and rj. 

Since ex is even and symmetric, r shifts to the left or to the right depending on 
its sign, and rj concentrates ex narrow or wide around the center of ex depending on 
its magnitude, the most interesting case for Ti is when the difference between x and 
r in ex is small. It is therefore sufficient to consider x G Z'^ and r G [n/2^n/2 + 1] 
where n is the dimension of Ti. Moreover, we can choose r/ > because is an even 
function. 

The following five pictures plot the point spectrum or eigenvalues denoted them as 
A of a finite Ti versus r G [n/2, n/2 + 1] with five different r] follows: 



3.5 h 

3h 
2.5- - 

2 - 
1.5 = 



0.5 

O' ' -u.o -U.J -U.J 

50 50.5 51 50 50.5 51 50 50.5 51 50 50.5 51 50 50.5 51 
11=0 r|=0.5 ri=1 r|=1.5 ii=2 

Figure 3. A vs r with n = 100 and r g [n/2, n/2 + 1] 
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The first graph 77 = shows that the range of the spectrum of Ti. is in [0, 4] which 
is similar as the discrete Laplacian operator where it is in [—4,0]. The second graph 
?7 = 0.5 shows that there is one isolated eigenvalue A = and the range of the other 
eigenvalues of 7i is smaller than it when rj = 0. The third graph where rj = 1 and 
the fourth graph where 77 = 1.5 show that there are two isolated eigenvalues and the 
range of all the other eigenvalues becomes smaller and smaller even than before. At 
last, the fifth graph shows that there are three isolated eigenvalues and the range of 
the eigenvalues besides those three eigenvalues of H is about [1.4,2.6]. However, at 
r = n/2 + 0.5 the third appearing eigenvalue seems to be still connecting with the 
spectrum, not totally isolated. Before I draw initial conclusions of critical values or 
intervals for r and 77, I like to show another five pictures that plot the eigenvalues A 
of Ti. versus 1/ cosh 17 where rj G [0, 2] with five different r as follows: 




0.5 10 0.5 10 0.5 10 0.5 10 0.5 1 

r=n/2 r=n/2+0.25 r=n/2+0.5 r=n/2+0.75 r=n/2+1 



Figure 4. A vs 1/ cosh 77 with n = 100 and rj e [0,2] 



The first graph where r = n/2 and the fifth graph where r = n/2 + 1 show that 
there are three isolated eigenvalues at 1/ coshrj ^ 0.3 and two isolated eigenvalues 
at 1/ coshr/ ^ 0.6. 0.3 corresponds approximately to rj = 2 and 0.6 is about 77 = 1. 
Moreover, the second graph where r = n/2 + 0.25 and the fourth graph where r = 
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n/2 + 0.75 show that there are still three isolated eigenvalues when 1/ cosh 77 f« 0.3 
and two isolated eigenvalues when 1/ cosh 77 f« 0.6, but now the third eigenvalue is 
not quite isolated from the set of other eigenvalues and seems to merge into the other 
sets of point spectrum. The third graph where r = n/2 + 0.5 show that there only 
exist two isolated eigenvalues. So, by Figure 3 and Figure 4 above, we can see that 
the third appearing eigenvalue does not seem to exist for all r e [n/2, n/2 + 1]. 

Therefore, we expect that the first and second isolated eigenvalues appear some- 
where at ?7 e (0, 1) and exist for all r e [^/2, n/2-|-l], and the third isolated eigenvalue 
appears somewhere at 77 e (1,2) and exist for for all r e [n/2,n/2 + 0.5) U (n/2 + 
0.5, n/2 + 1]. By approximating infinite H. from finite these three eigenvalues will 
remain in the point spectrum and all other eigenvalues will belong in the continuous 
spectrum. 

You may ask that how numerical calculations for finite size systems can give insight 
in the infinite size systems. Since Ti. is defined on ^^(Z), we can focus on a subspace of 
£^(Z) such that H is finite dimensional. Eigenvalues and eigenvectors of finite H can 
be calculated by numerics and we can increase the size of Ti since the eigenvectors 
are defined on £^ (Z) . The spectrum of infinite dimensional H. can then be predicted 
by observing it of finite dimensional H.. 

Moreover, we know the structure of the spectrum of the discrete Laplacian operator 
and is a perturbed case of it, we can study the spectrum of the discrete Laplacian 
operator in order to understand the spectrum of H. based on a theorem called the 
compact perturbation theorem below. 

Theorem 2. Let T e -B(X) and let AheT — compact. Then T and T + A have the 

same essential spectrum. 

This theorem is rewritten from 'Perturbation Theory for Linear Operators' by Kato 
[K]. 

Ti. is constructed basically by reversing the sign of the discrete Laplacian operator 
and adding a compact perturbation term e^. The operator e^I in Ti. is in fact a 



16 

compact perturbation. To prove this statement, we can first define T — [e^ — 2)X 
because we intend to show that e^l is perturbed from 2X which is the diagonal of 

the discrete Laplacian operator with inverse sign. We can then define an operator 7^ 
with finite-dimensional range such that 7^ is compact [HN] and 
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If lim^^oo ll'^t ^ '^11 = 0) then T is compact where || ■ || denotes the operator norm. 
Now, we let v e ^^(Z) and 



< max|6.-2|(5^z.,2)V2 

|a;|>n 

= max|ea; - 2|||i/||£2(z) 

|a;|>n 

So, lim„^oo ||7^ — T|| = because lim3.^oo — 2| = is as shown in Figure 1. 

Thus, e^-X in 7i is a compact perturbation and Ti and the discrete Laplacian op- 
erator have the same essential spectrum. The spectrum of the discrete Laplacian 
operator is purely continuous and is contained in the interval [—4, 0] and the point 
spectrum of finite at = is contained in [0, 4] based on my numerics. When 
77 = 0, 7i is different from the discrete Laplacian operator by a sign. Moreover, my 
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numerical results presented in the next section shows that at 77 = the maxima of the 
square normalized components of all the eigenvectors of the three isolated eigenval- 
ues decrease as the matrix size n oiTi increase and we expect that the eigenvectors 
converge to zero for infinite n. So, there are no eigenvalues when 77 = for infinite 
dimensional Ti. and they all converge to the set of the continuous spectrum. This fact 
agrees with the case of the discrete Laplacian operator. Therefore, when we change r 
and 77 of a finite n-dimensional H, we expect that the three isolated eigenvalues will 
converge to the set of the point spectrum and all other eigenvalues will converge to 
the set of point spectrum as n goes to infinity. Finally, the residual spectrum of the 
discrete Laplacian operator and Ti are empty because they are both bounded and 
self-adjoint [HN]. 

The mathematical argument that determines the continuous spectrum in the infi- 
nite chain limit of the discrete Laplacian operator is presented below. 

Let T : L2(T) e'^(R) denote the Fourier transform and J^* be the inverse Fourier 
transform. T is an unitary operator. For proofs, see Hunter and Nachtergaele [HN]. 
Moreover, consider only now that A = S + S* — 2X is the discrete Laplacian operator. 

It follows that T*AT — T> where P is a diagonal matrix that consists of spectral 
elements of A on the diagonal since A is self-adjoint. Moreover, if f{x) e L^(T) 
where T = [0, 27r] and ^ e then 
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T*/^Tf{x) = T*{S + S* -2X)J^f{x) 

— ^*{fn-l + fn+1 — '^fn) 

neZ neZ ngZ 

(2)^E/»^- 

= (e^^ + e-*^ - 2)/(x) 
= (2cos(a;) -2)/(a;) 

= ^^/(^) 

It follows that (7(A) = a{V) = e [0,27r]} = [-4,0] and a{A) is entirely 

continuous. 

Moreover, the bottom of the continuous spectrum of Ti. is given by a formula 
2(1 — A~^) [MN] which agrees with my numerics shown in Figure 3 above. The 
smallest eigenvalue, or the first isolated eigenvalue of H in our case, is also called 
the "ground state energy" and its corresponding eigenfunction is called the "ground 
state" which is the configuration of the system with the smallest total energy, and 
the other eigenvalues correspond to discrete excitations of the system [D] . 

As proved by Michoel and Nachtergaele [MN] , the eigenfunction of the first isolated 
eigenvalue of H is V^^^ = — ^ e f(Z). 
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The fact that this eigenfunction belongs to £^(Z) can proved as follows: 



^ ' cosh(77(a; - r)) ' ^ e2^(^-'') + e-2'?(^-'-) + 2 
and 

4 46^'"^ 



< 



g277(a;— r) ^ 2r;(a;— r) _|_ 2 — g2f?(x— r) g2?73; 

Since ^^.g^ is a geometric series, by comparison test, we can conclude that 



E 



< oo 



cosh (77(3; — r)) 

and so K^'^ e 

Moreover, the fact that V^^ is indeed an eigenvector of JiV^^ — for the zero 
mode can be verified as follows: 



cosh (77(0; — r))' 



^^^cosh(r7(x — r))' 
1 . 1 



+ 



-) 



A cosh(?7(a; — 1 — r)) cosh(?7(x + 1 — r)) ' 

2 cosh(?7(a; — r))^ 1 



-)(- 



'cosh(77(x — 1 — r)) cosh(77(a; + 1 — r)) cosh(r7(a; — r)) 
1.1 1 



) 



cosh(r7 



cosh(?7 
1 



cosh(r7 



cosh(?7 
1 



cosh(?7 



cosh(r7 



cosh(?7 




+ 



'cosh(?7(a; — 1 — r)) cosh(r7(a; + 1 — r)) 

2 cosh(?7(a; — r)) 



) 



X — 1 — r)) cosh(?7(a; + 1 — r)) 
cosh(r7(a; + 1 — r)) + cosh(?7(,T — 1 — r)) 
cosh(r7(a; — 1 — r)) cosh(?7(a; + 1 — r)) 
2 cosh(77(a; — r)) 



a; — 1 — r)) cosh(?7(a; + 1 — r)) 
^ ((2cosh(?7))(2cosh(a; -r)))/2 
cosh(?7(a; — 1 — r)) cosh(?7(a; + 1 — r)) ' 
2 cosh(?7(x — r)) 



a; — 1 — r)) cosh(r7(a; + 1 — r)) 
2 cosh(r7(a; — r)) 



a; — 1 — r)) cosh(r7(a; + 1 — r)) 



20 

The following section explains my MATLAB programs, and it also shows my nu- 
merical results for all the three isolated eigenvalues in order. 
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4. The Eigenvalue Problem 

There exist two isolated eigenvalues of H. when 77 > and r G [n/2,n/2 + 1]. The 
third isolated eigenvalue appears when 77 > 1 and the interval of r where it does not 
exist becomes smaller and smaller when 77 becomes larger and larger. Moreover, the 
third eigenvalue does not exist at r = n/2 + 0.5. These numerical results are based 
on three programs which are explained in the next paragraph. 

The first program named "spectrum.m" plots the entire set of eigenvalues of an 
n-dimensional Ti, versus r e [n/2,n/2 -\- 1\ and -^^^ where 77 e [0,2] separately as 
shown in Figure 3 and Figure 4 above. This program gives the picture of the lo- 
cations in which there exist isolated eigenvalues. With various r and 7] in 7Y, the 
second program named "mainresults.m" graphs the eigenvectors of all the three iso- 
lated eigenvalues versus the components of the eigenvectors, and this program also 
calculates the maxima of the square normalized components in the eigenvectors of 
those eigenvalues. The first part of "mainrestuls.m" shows the distributions of the 
eigenvectors at different r and 77. The second part compares the maxima described 
above to see when they do not decrease for certain r and 77 as n increases which 
imply the eigenvectors are nonzero and the existence of eigenvalues. The last pro- 
gram named "dccreasingratcs.m" is designed only for the third appearing eigcnvahic 
because it does not exist for all r e [7i/2, 7i/2 -|- 1]. I claim that the third eigenvalue 
does not exist when r = 7^/2 -|- 0.5 and so I investigate the decreasing rate of three 
selected components of the eigenvector of the third eigenvalue. We let denote 
an n-dimensional eigenvector at components x = 77-/4, 77,/2, 3rz/4 and assume that 
< because the square normalized component is decreasing as shown on 
Figure 7 below, then defined to be the decreasing rate. We are 

expecting to see that a— >Oasr— >(n/2-|- 0.5). 

The codes of all the programs are described in Appendix B below. The follow- 
ing three subsections show some selected numerical results for the three eigenvalues 
denoted them as Ai, A2, and A3, respectively. 
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4.1. The first isolated eigenvalue. The pictures below graph the eigenvector de- 
noted as vi^^ with respect to Ai versus the components denoted the index as x of 
Vx^'^ when n = 100. From top to bottom, r increases from n/2 to n/2 + 1 by 0.25 
each time, and from left to right, rj increases from to 2 by 0.5 each time. 




50 1 00 



50 1 00 



50 1 00 



50 1 00 



50 100 



50 100 



50 100 



50 100 



50 1 00 



50 100 



50 1 00 



50 1 00 



50 100 50 100 50 100 50 100 50 100 
ri=0 Ti=0.5 11=1 n=1.5 n=2 



Figure 5. vP vs X with n = 100 



As you can see from Figure 5 above, when rj > 0, Vx^^ starts to concentrate around 
X = n/2 for all r G [n/2, n/2 + 1] and there exists an absolute maximum at x = n/2 on 
all the graphs. This maximum shows that Vx'^^ is nonzero although most components 
of vP are zero and implies the existence of the eigenvalue Ai for infinite 7i. Similarly 
for the cases of larger n which are not shown here, Vx^^ has an absolute maximum at 
X = n/2 and is zero almost everywhere else. 

Next, 1 want to justify that if the absolute maximum mentioned above does not 
decrease for all rj > and r G [n/2, n/2 + 1] as n increases. 

The table below shows the maxima of the square normalized components in vj^^ 
at n = 200 and rj G [0.2] where t] is increased by (2 — 0)/4 = 0.5 every time, and r 
ranges from n/2 to n/2 + 1 and is increased by 0.25 each time. 
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77 = 


T] = 0.5 


1] = 1 


1] = 1.5 


r] = 2 


r = n/2 


0.009950 


0.250000 


0.498981 


0.723493 


0.874100 


r = n/2 + 0.25 


0.009950 


0.246134 


0.470008 


0.653751 


0.788057 


r = n/2 + 0.5 


0.009950 


0.235004 


0.394028 


0.464358 


0.488313 


r = n/2 + 0.75 


0.009950 


0.246134 


0.470008 


0.653751 


0.788057 


r = n/2 + 1 


0.009950 


0.250000 


0.498981 


0.723493 


0.874100 



The table below shows these maxima at n = 300 and the rest is the same as above. 





rj = 


ri = 0.5 


1] = 1 


1] = 1.5 


r] = 2 


r — n/2 


0.006644 


0.250000 


0.498981 


0.723493 


0.874100 


n/2 + 0.25 


0.006644 


0.246134 


0.470008 


0.653751 


0.788057 


n/2 + 0.5 


0.006644 


0.235004 


0.394028 


0.464358 


0.488313 


r = n/2 + 0.75 


0.006644 


0.246134 


0.470008 


0.653751 


0.788057 


r = n/2 + 1 


0.006644 


0.250000 


0.498981 


0.723493 


0.874100 



By comparing the maxima of the square normahzed components in V-} on these 
two tables at the same values of r and r} with two different n, it shows that the 
maxima remain unchanged besides when 77 = 0. When rj — 0, the maxima decrease 
as n increases for all five different r. So the maxima begin to be nondecreasing 
somewhere at rj E (0, 0.5). My next step is to compare the maxima between n = 300 
and n = 400 with 77 e [0, 0.5] and increase 77 by (0.5 — 0)/4 — 0.125 each time and keep 
the same r as above. The following table shows the procedure of how I narrow down 
the interval of 77 with different ni, n2 where ni < n2 and 771, 772 where 77 e [771, 772]. 
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ni 


n2 


Vi 


V2 


200 


300 





2 


300 


400 





0.5 


400 


500 





0.125 


500 


600 





0.0625 


600 


700 





0.03125 



For the tables from n — 700 to n — 1000 which are not shown here when rj e 
[0, 0.03125], the interval of rj starting at zero becomes smaller and smaller beginning 
with r = n/2 and r = n/2 + 1 where the maxima of the square normalized components 
in Vx^^ do not decrease, followed by r = n/2 + 0.5, and then r — n/2 + 0.25 and 
r = n/2 + 0.75. Since the spectrum for infinite H is purely continuous when 77 = 
because it is only different from the discrete Laplacian operator by a sign and the 
point spectrum becomes nonempty for an rj in (0, 0.3125) which is an open set, the 
critical value for r] where Ai exists for infinite 7i is believed to be at 77 = 0. 

4.2. The Second isolated eigenvalue. A2 exists in a very similar fashion as Ai. 
The pictures below draw the eigenvector denoted as Vx with respect to A2 versus 
the components denoted the index as x of Vx'^^ when n = 100. Prom top to bottom, 
r is increased from n/2 to n/2 + 1 by 0.25 each time, and from left to right, 77 is 
increased from to 2 by 0.5 each time. 

As appeared on Figure 6 below, Vx"^^ starts to be concentrated obviously around 
X — n/2 like Vx^^ when 77 > 0.5 for all r e [n/2, n/2 + 1] and there also exists an 

absolute maximum at x = 7i/2 on those graphs. By graphing the same pictures with 

(2) 

larger n, there always exists an absolute maximum in . Thus, the eigenvector will 
remain nonzero and A2 will not disappear for large n. 

The two tables presented below show the maxima of the square normalized compo- 
nents in yi^^ at 71 = 200 and n — 300, and 77 e [0.2] where 77 increases by (2— 0)/4 = 0.5 
every time, and r ranges from 7i/2to7i/2-|-l and increases by 0.25 each time. 
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Figure 6. vs X with n = 100 





= 


T] = 0.5 


f] = 1 


T] = 1.5 


7] = 2 


r = n/2 


0.009950 


0.032095 


0.159968 


0.375804 


0.459508 


r = n/2 + 0.25 


0.009950 


0.032057 


0.171671 


0.526199 


0.772677 


r = n/2 + 0.5 


0.009950 


0.031823 


0.172972 


0.368510 


0.476093 


r = n/2 + 0.75 


0.009950 


0.032057 


0.171671 


0.526199 


0.772677 


r = n/2 + 1 


0.009950 


0.032095 


0.159968 


0.375804 


0.459508 





r] = 


7] = 0.5 


f] = 1 


T] = 1.5 


7] = 2 


r = n/2 


0.006644 


0.031971 


0.159968 


0.375804 


0.459508 


r = n/2 + 0.25 


0.006644 


0.031933 


0.171671 


0.526199 


0.772677 


r = n/2 + 0.5 


0.006644 


0.031699 


0.172972 


0.368510 


0.476093 


r = n/2 + 0.75 


0.006644 


0.031933 


0.171671 


0.526199 


0.772677 


r = n/2 + 1 


0.006644 


0.031971 


0.159968 


0.375804 


0.459508 



The maxima shown on the tables above appeare to be nondecreasing when 77 > 0.5 
and my next step is to focus on 7] E [0, 1] and compare the maxima of the square 
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normalized components of Vx"^^ between n — 300 and n — 400. My complete process 
is shown in the table below: 





n2 


Vi 


V2 


200 


300 





2 


300 


400 





1 


400 


500 





0.75 


500 


600 





0.5625 



When n = 600,700,800,900,1000 for r] e [0,0.5625], same as in the previous 
situation but much slower, the interval of r] starting at zero where the maxima do 
not decrease becomes smaller and smaller beginning with r = n/2 and r = n/2 + 1, 
followed by r = n/2 + 0.5, and then r = n/2 + 0.25 and r = n/2 + 0.75. Although the 
convergence is slower in this case, the critical value of rj where A2 exists for infinite 
H is still at = 0. 

4.3. The third isolated eigenvalue. Finally, A3 exists in a situation that is quite 
different from it of Ai and A2. The pictures shown below graph the eigenvector 
denoted as Vx^^ with respect to A3 versus the components denoted the index as x of 
Vx^^ when n — 100. Prom top to bottom, r is increased from n/2 to n/2 + 1 by 0.25 
each time, and from left to right, 77 is increased from to 2 by 0.5 each time. 

As shown on Figure 7 below, for r = n/2 and r = n/2 + 1, V^^^'' do not seem to have 
their concentrated regions or absolute maxima until rj > 1, and for r = n/2 + 0.25 
and r = n/2 + 0.75, the absolute maxima appear when rj > 1.5. At r = n/2 + 0.5, 
Vx^^ do not focus around x = n/2 when n = 100, and the maxima become smaller 
and smaller when n increases. Thus, A3 may not be an eigenvalue when r = n/2 + 0.5 
because Vx^^ may be zero at r — n/2 + 0.5 for infinite H. 

The two tables presented below show the maxima of the square normalized com- 
ponents in Vx^^ at n = 200 and n = 300, and rj e [0.2] where rj is increased by 
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r=n/2+ 
0.25 



r=n/2+ 
0.5 



r=n/2+ 
0.75 




50 100 50 100 50 100 50 100 50 100 
71=0 11=0.5 11=1 n=1-5 11=2 



Figure 7. vs X with n=100 



(2 — 0)/4 = 0.5 every time, and r ranges from ?2/2 to n/2 + 1 and is increased by 0.1 
each time. 





ri = 


r] = 0.5 


rj = 1 


r] = 1.5 


7] = 2 




r = n/2 


0.012558 


0.012896 


0.119863 


0.306798 


0.389453 


r 


= n/2 + 0.1 


0.011965 


0.011392 


0.116918 


0.548307 


0.767640 


r 


= n/2 + 0.2 


0.011402 


0.010295 


0.014707 


0.476112 


0.706397 


r 


= n/2 + 0.3 


0.010892 


0.010118 


0.015720 


0.195049 


0.508429 


r 


= n/2 + 0.4 


0.010433 


0.010182 


0.013939 


0.018121 


0.042290 


r 


= n/2 + 0.5 


0.010005 


0.009957 


0.009903 


0.009841 


0.009763 


r 


= n/2 + 0.6 


0.010433 


0.010182 


0.013939 


0.018121 


0.042290 


r 


= n/2 + 0.7 


0.010892 


0.010118 


0.015720 


0.195049 


0.508429 


r 


= n/2 + 0.8 


0.011402 


0.010295 


0.014707 


0.476112 


0.706397 


r 


= n/2 + 0.9 


0.011965 


0.011392 


0.116918 


0.548307 


0.767640 


r = n/2 + 1 


0.012558 


0.012896 


0.119863 


0.306798 


0.389453 
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r] = 


ri = 0.5 


1] = 1 


1] = 1.5 


r] = 2 


r 


= n/2 


0.008375 


0.008632 


0.119863 


0.306798 


0.389453 


r — 


n/2 + 0.1 


0.007978 


0.007622 


0.116918 


0.548307 


0.767640 


r — 


n/2 + 0.2 


0.007602 


0.006883 


0.010120 


0.476112 


0.706397 


r = 


n/2 + 0.3 


0.007261 


0.006757 


0.010572 


0.195049 


0.508429 


r = 


n/2 + 0.4 


0.006955 


0.006798 


0.009331 


0.012241 


0.039502 


r — 


n/2 + 0.5 


0.006669 


0.006648 


0.006624 


0.006596 


0.006560 


r — 


n/2 + 0.6 


0.006955 


0.006798 


0.009331 


0.012241 


0.039502 


r — 


n/2 + 0.7 


0.007261 


0.006757 


0.010572 


0.195049 


0.508429 


r = 


n/2 + 0.8 


0.007602 


0.006883 


0.010120 


0.476112 


0.706397 


r — 


n/2 + 0.9 


0.007978 


0.007622 


0.116918 


0.548307 


0.767640 


r - 


= n/2 + l 


0.008375 


0.008632 


0.119863 


0.306798 


0.389453 



As presented on the two tables above, at r = n/2, r = n/2 + 0.1, r = n/2 + 0.9, 
and r = n/2 + 1, the maxima of the square normalized components in Vx^^ do not 
decrease when rj > 1. Moreover, at r = n/2 + 0.2, r = n/2 + 0.3, r = n/2 + 0.7, 
and r = n/2 + 0.8, the maxima do not decrease when rj > 1.5. Lastly, the maxima 
decrease at r = n/2 + 0.4, r = n/2 + 0.5, and r = n/2 + 0.6 for all rj e [0, 2]. Based 
on these two tables, we may expect that A3 exists when rj > 1 while the interval of r 
where A3 does not exist is becoming smaller and smaller as 1] is becoming larger and 
larger. However, up to now there is not enough argument to support the claim that 
A3 does not exist at r = n/2 + 0.5. The pictures below test this claim. 

Figure 8 below shows the decreasing rates mentioned above at x — n/4, x — n/2, 
and a; = 3n/4 of Vx^^ with n = 100, n = 200, and n = 300. As you can see that the 
decreasing rate a approaches zero as r tends to n/2 + 0.5. When the decreasing rate 
is zero, it implies that the eigenvector at that component flattens out which does not 
increase or decrease. Therefore, Vx^^ is expected be zero at r = n/2 + 0.5 when n 
goes to infinity. 
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0.4 p 
0.3 
0.2 
0.1 



50 

0.8 p- 

0.6 

0.4 

0.2 

— 
50 



50.5 
n=100 



100 

3| — 



100 

0.8 I 

0.6 
0.4 
0.2 



100.5 
n=200 



101 150 

6| 



101 150 

3| 



150.5 
n=300 



Figure 8. a vs r with rj = 2 

The graphs of the eigenvectors gives a conclusion that A3 appears when rj > 1 
while the tables of maxima shows that it happens when t] > 1. So from my numerical 
results, I conclude that the critical value of rj where A3 exists at the borders of 
r G [n/2,n/2 + 1] is rj = 1, and the interval of r where A3 does not exist becomes 
smaller and smaller as 77 becomes larger and larger and does not exist atr = n/2 + 0.5. 
The critical value of rj means that when 77 < 1, A3 does not exist for any r, but A3 
starts to exist when 77 > 1 for certain interval of r. 
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5. Conclusion 

My numerical results predict that the spectrum of infinite dimensional Ti contains 
not only the continuous spectrum but also the point spectrum. The residual spectrum 
is empty because 7i is a bounded and self-adjoint operator. The point spectrum 
of infinite dimensional H. consists of three eigenvalues with particular regions of r 
and 7]. The results come from using finite dimensional Ti to approximate infinite 
dimensional Ti. The reason why we can use finite H. to estimate infinite H. is based 
on the compact perturbation theorem because we are already aware of the case of the 
discrete Laplacian operator which is also an unperturbed case of H.. Since for finite 
dimensional 7i the spectrum contains only the point spectrum or a set of eigenvalues, 
I focus on the three eigenvalues that are isolated from the set of the other eigenvalues 
because those three isolated eigenvalues are expect to remain in the point spectrum 
for infinite dimensional Ti.. These eigenvalues still exist for infinite H. when their 
corresponding eigenvectors are nonzero depending on r and Tj. It is sufficient to start 
with re [n/2,n/2-|-l] and > because of the properties of the diagonal term 
in Ti.. Based on the numerics produced by my programs, I conclude that the first 
(Ai) and second (A2) eigenvalues exist when 77 > and for all r e [n/l, n/2 + 1]. The 
third cigcnvahic (A3) exists when 77 > 1 and the interval of r where A3 does not exist 
becomes smaller and smaller as r] becomes larger and larger. Moreover, A3 does not 
exist at r = n/2 + 0.5. 

Some further discussions are listed as follows: First, larger matrix size n and pa- 
rameter T] are needed for a more complete solution of the existence of A3 because 
my programs can only conduct the calculations up to n = 300 and rj = 2. Second, 
Michoel and Nachtergaele [MN] have proved that v}^^ — cosh{rl^{x-r)) ■ ^'-'^ ^^"^^ 
and Vx^\ eigenfunctions hke Vx^^ may or may not exist. Third, a detailed analytic 
proof for using finite dimensional TC to approximate infinite dimensional TC is needed 
because in my paper I state that this path is possible because of our knowledge of 
the discrete Laplacian operator. All these questions are very interested to be studied 
in future research. 
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Appendix A. Definitions 



A.l. The definition of a Hilbert space. 



Definition 3. An inner product on a complex linear space 'K is a map 



(•,•) :XxX^C 



such that, for all x,y,z E'K. and X, Ijl E C; 

(1) {x, \y + fiz) = X{x, y) + z) (linear in the second argument); 

(2) {y-ix) — {x,y) (Hermitian symmetric); 

(3) {x,x) > (nonnegative) ; 

(4) {x,x) = if and only if x — (positive definite); 

We call a linear space with an inner product an inner product space or a pre-Hilbert 



Definition 4. A Hilbert space is a complete inner product space. 
A.2. The definition of £p(Z) and Lp{W). 

Definition 5. For 1 < p < oo, the sequence space £^(Z) consists of all infinite 
sequences x — {xn)'^=-(x> ■^^'^^ ^^^^ 



Definition 6. Suppose that I < p < oo. We denote by U\M.) the set of Lebesgue 
measurable functions f -.M. —>■ R. (or C) such that 



where the integral is a Lebesgue integral, and we identify functions that differ on a 
set of measure zero. 



space. 



oo 
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A. 3. The definition of linear, bounded, compact, and self-adjoint operators. 

Definition 7. A linear map or linear operator T between real (or complex) linear 
spaces X, Y is a function T : X — > Y such that 

T{Xx + fxy) ^ XTx + fj,Ty 

for all X , II & R (or C) and x, y & X.. A linear map T : X ^ Y is called a linear 
transformation o/X, or a linear operator on X. 

Definition 8. Let X and Y be two normed linear spaces. We denote both the X and 
Y norms by \\ ■ \\. A linear map T : X ^ Y is bounded if there is a constant M > 

such that 

\\rx\\ < M\\x\\ 

for all X & 'K. If no such constant exists, then we say that T is unbounded. 

Definition 9. A linear operator T : X ^ Y is compact if and only if every bounded 
sequence (a;„) in X has a subsequence (a;„^) such that {Txnf.) converges in Y. 

Definition 10. A bounded linear operator A : H ^ H on a Hilbert space H is 
self-adjoint if and only if 

{x, Ay) = {Ax, y) 

for all y e H 

A. 4. The definition of kernel and range of an operator. 

Definition 11. Let T : X ^ Y &e a linear map between linear spaces X, Y. The 
null space or kernel of T , denoted by kerT , is the subset of X defined by 

kerT = {x & X|Ta; = 0}. 

The range of T , denoted by ranT, is the subset of Y defined by 

ranT — {y & Y\3x e X|Tx = y}. 
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A. 5. The definition of the Fourier transform and inverse Fourier transform. 

Definition 12. The periodic Fourier transform T : L^(T) £^(Z) that maps a 
function to its sequence of Fourier coefficients is defined by 



n=—oo 



such that 

fix) = - 

/27r 

where T= [0,27r]. 

Definition 13. The inverse Fourier transform T* : £^(Z) — >• L^(T) is defined by 



^ oo 



n=—oo 



A. 6. The definition of unitary operators. 

Definition 14. A linear map U : Hi — > H2 between real or complex Hilbert spaces 

Hi and H2 is said to he orthogonal or unitary, respectively, if it is invertible and if 

for all x,y & Hi. Two Hilbert spaces Hi and H2 are isomorphic as Hilbert spaces if 
there is a unitary linear map between them. 
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Appendix B. Programs 

The codes of my three programs: "spectrum.m" , "maimesults.m", and "decreasin- 
grates.m" are presented below, "spectrum.m" is in the first box, "mainresutls.m" is 
in the second box and third box, and "decreasingrates.m" is in the fourth box. The 
percentage signs appeared in the codes represent comment for descriptions and not 
functioning for certain tasks because "spectgrum.m" plots A versus r and rj separately 
and "mainresults.m" graphs Vx^\ Vx'^\ Vx'^^ and calculates the maxima of the square 
normalized components of those eigenvectors at different time as well. 
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%plotting the eigenvalues 

x=100; A=zeros(x,x); 
k=l ; %subplot index 

%for eta=0:0.5:2 %lambda vs r 
for r=x/2:0.25:x/2+l %lambda vs eta 
%for r=x/2:0.0 1 :x/2+l %for lambda vs r 
for eta=0:0. 05:2 %lambda vs eta 
%calculating v 
for y=l :x 

epsilon(y)=2*(cosli(eta*(y-r)))^2/(cosh(eta*(y-l-r))*cosh(eta*(y+l-r))); 
end 

%assignning entries 
for i=l ;x 

A(i,i)=epsilon(i); 
end 

for i=l:x-l 

A(i,i+l)=-l/cosh(eta); 

end 

for i=2;x 

A(i,i- 1 )=- 1 /cosh(eta) ; 
end 

%diagnalizing A 
V=eig(A); 

% %plotting lambda vs r 

% subplot(l,5,k) 

% plot(r,V) 

% hold on 

%plotting lambda eta 
subplot(l,5,k) 
plot(l./cosh(eta),V) 
hold on 

end 

k=k+l; 

%r=x/2; %for lambda vs r 
eta=0; % lambda vs eta 
end 



Figure 9. spectrum.m 
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%verifying the three isolated eigenvalues 

%matrix size 
x=100; 

j=l ; %subplot index 
k= 1 ; %r index 
1=1; %eta index 

%interval of eta 
a=0; b=2; 

% defining A 
A=zeros(x,x); 

for r=x/2:0.25:x/2+l 
for eta=a:(b-a)/4:b 

%calculating epsilon 
y=[l:l:x]'; 

epsilon(y)=2*(cosh(eta.*(y-r))).^2./(cosh(eta.*(y-l-r)).*cosh(eta.*(y+l-r))); 

%assignning entries 
for i=l :x 

A(i,i)=epsilon(i); 
end 

for i=l :x-l 

A(i,i+l)=-l/cosh(eta); 
end 

for i=2:x 

A(i,i- 1 )=- 1 /cosh(eta) ; 
end 

%finding eigenvalues and eigenvectors of A 
[V,D]=eig(A); 

%plotting a normalized eigenvector vs its component 
xx=l :1 :x; 

yy=V(:,3)/norm(V(:,3)); 
subplot(5,5,j) 
plot(xx,yy) 
hold on 

j=j+i; 

% %calculating the max 

% M(k,l) = max((V(:,2)/norm(V(:,2))).'^2); 

% 1=1+1; 



Figure 10. mainresults.m 
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end 

1=1; 

k=k+l; 
end 



% [x a b] %disp laying the matrix size and the interval of eta 

% 

% %displaying the max 
% for c=l:k-l 

% disp(sprintf('%4f & %4f & %4f & %4f & %4f \V,M(c,l),M(c,2),M(c,3),M(c,4),M(c,5))) 
% end 



Figure 11. mainresults.m cont. 
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%observmg the tail part of the eigenvector for the three isolated eigenvalues 

x=100; %initial matrix size 
1=1 ; %subplot index 
eta=2; 

for k=x/4:x/4:3*x/4 
for x=100:100:300 
A=zeros(x,x); 
for r=x/2:0.02:x/2+l 

%calculating epsilon 
y=[l:l:x]'; 

epsilon(y)=2*(cosh(eta.*(y-r))).^2./(cosh(eta.*(y-l-r)).*cosh(eta.*(y+l-r))); 

%assignning entries 
for i=l :x 

A(i,i)=epsilon(i); 
end 

for i=l :x-l 

A(i,i+l)=-l/cosh(eta); 
end 

for i=2:x 

A(i,i-l)=-l/cosh(eta); 
end 

%diagnalizing A 
[V,D]=eig(A); 

%defining the decreasing rate a 
a=-log((V(k,3 )/norm(V(: ,3 )))^2)/k; 

%plotting r vs a 

subplot(3,3,l) 
plot(r,a) 
hold on 
end 
1=1+1; 
end 
end 



Figure 12. decreasingrates.m 
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